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SUMMARY 


Diffusion calculations were performed to establish the conditions under 
which concentration dependence of the diffusion coefficient was important in 
single-, two-, and three-phase binary alloy systems. Finite-difference solu- 
tions for each type of system were obtained using diffusion coefficient varia- 
tions typical of those observed in real alloy systems. Solutions were also 
obtained using average diffusion coefficients determined by taking a logarithmic 
average of each diffusion coefficient variation considered. The solutions for 
constant diffusion coefficients were used as references in assessing the effects 
of diffusion coefficient variations. Calculations were performed for planar, 
cylindrical, and spherical geometries in order to compare the effect of diffu- 
sion coefficient variations with the effect of interface geometries. 

Diffusion coefficient variations in single-phase systems and in the major- 
alloy phase of two-phase systems were found to effect the kinetics of diffusion 
as strongly as the interfacial geometry of the diffusion couple. Concentration 
dependence of the diffusion coefficient in the minor-alloy phase of a two-phase 
system was found to have only an initial transient effect on the diffusion 
kinetics. In three-phase systems, the intermediate phase did not increase in 
thickness if the diffusion coefficient of the intermediate phase was smaller 
than the diffusion coefficient of the major-alloy phase. Under these condi- 
tions, the three-phase diffusion problem could be treated as a two-phase problem. 
However, if the diffusion coefficient of the intermediate phase was larger than 
the diffusion coefficient of the major-alloy phase, the intermediate phase grew 
rapidly at the expense of the major and minor phases. In most of the cases con- 
sidered, the diffusion coefficient of the major-alloy phase was the key param- 
eter that controlled the kinetics of interdiffusion. 


INTRODUCTION 

Solutions to the diffusion equation for a variety of initial and boundary 
conditions have been reported in the literature (refs. 1 to 7) . Most of these, 
however, are restricted to planar geometry, are applicable only for infinite or 
semi-infinite systems, or require that the diffusion coefficient be constant. 

In treating diffusion in binary alloy systems, the assumption that the diffu- 
sion coefficient is constant is not always valid. Diffusion coefficient data 
for several single-phase and multiphase binary alloy systems show that an order 
of magnitude variation in the diffusion coefficient within a phase is not 
uncommon and in some cases, the diffusion coefficient may vary by as much as 
two orders of magnitude (refs. 8 to 12) . A determination of the relative impor- 
tance of the types and magnitudes of the diffusion coefficient variations found 
in real systems would provide an indication of the magnitude of the error asso- 
ciated with treating the diffusion coefficient as a constant. 

The present study was undertaken to determine the relative importance of 
diffusion coefficient variations in single-, two-, and three-phase binary alloy 



systems with planar, cylindrical, or spherical interfaces. Numerical solutions 
to the governing diffusion equations were obtained for various concentration- 
dependent diffusion coefficients as well as for two average (constant) diffusion 
coefficients. Two different methods of averaging were considered to determine 
the averaging method best suited for different types of diffusion coefficient 
variations. Solutions were also obtained for planar, cylindrical, and spherical 
geometries to compare the effects of interface geometries with those caused by 
concentration-dependent diffusion coefficients. 
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SYMBOLS 

elements which constitute binary alloy system 

concentration of sample, atomic fraction of B 

average concentration of sample, atomic fraction of B 

initial • concentration of 3-phase, assumed to be unity in this study, 
atomic fraction of B 

initial concentration of a-phase, assumed to be zero in this study, 
atomic fraction of B 

chemical diffusion coefficient (concentration dependent) , m /sec 
chemical diffusion coefficient of pure A, m 2 /sec 
chemical diffusion coefficient of pure B, nr /sec 

designates a type of log linear D variation, D increases going 
away from interface, m^/sec 

designates a type of log linear D variation, D decreases going 
away from interface, m 2 /sec 

designates a concentration-independent D whose value equals the 
logarithmic average of D f or D + variations considered in this 
study, m 2 /sec 

designates a concentration-independent D whose value equals the 
integral average of D + or D + variations considered in this 
study, m 2 /sec 

amount of B in 3-phase remaining at time t normalized to initial 
value for two- and three-phase systems; or amount of B remaining 
in original B-rich region at time t normalized to initial value 
for single-phase system 

diffusion flux, atoms/rn 2 /sec 
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L total thickness or diameter of specimen, m 

£ initial thickness or diameter of B-rich region, m 

R normalized distance, r/(L/2) 

r distance from center of B-rich region, m 

s interfacial surface area of B-rich material, m 2 


t 


V 


3 


V 


a, 6 


Y 


6 


3a 


0D 

K/2 

Cl/2 

C 2 /2 


time, sec 

volume fraction of B-phase in specimen at equilibrium; other sub- 
scripts have similar meaning 

volume of B-rich material, m 2 

designations for terminal phases rich in A and B, respectively 

designation for intermediate phase 

difference between solubility limits of 3- and a-phases, C^ a - C a ^ 
(other double subscripts have similar meaning) , atomic fraction 
of B 

temperature of interest, K 

distance from center of B-phase to aB interface in two-phase system, 
positive superscript designates a side of interface and negative 
superscript designates 3 side of interface, m 

distance from center of 3-phase to BY interface in three-phase 
system, positive superscript designates y side of interface and 
negative superscript designates 3 side of interface, m 

distance from center of B-phase to ay interface in three-phase 
system, positive superscript designates a side of interface and 
negative superscript designates y side of interface, m 


Superscripts: 
d A-rich phase 


B B-rich phase 

y intermediate phase 

* denotes equilibrium value 
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Subscripts : 


aB 

designates 

location 

in 

a-phase 

at 

aB 

interface 

of 

two-phase system 

ay 

designates 

location 

in 

a-phase 

at 

ay 

interface 

of 

three-phase 

system 

Ba 

designates 

location 

in 

B-phase 

at 

aB 

interface 

of 

two-phase system 

By 

designates 

location 

in 

B-phase 

at 

By 

interface 

of 

three-phase 

system 

ya 

designates 

location 

in 

y-phase 

at 

ay 

interface 

of 

three-phase 

system 

yB 

designates 

location 

in 

y-phase 

at 

By 

interface 

of 

three-phase 

system 


MATHEMATICAL ANALYSIS 

The relationships between the binary phase diagrams and concentration pro- 
files produced by diffusion between pure A and pure B after annealing for 
some intermediate time at temperature 0 D in single-, two-, and three-phase 
binary alloy systems are shown in figure 1. Figure 1(a) shows a continuous 
diffusion profile for a binary system A-B with complete solid solubility 
between A and B. Diffusion in a two-phase system with limited solid solu- 
bility produces a concentration profile (fig. 1(b)) with discontinuity at the 
oc$ interface equal to the difference between the solubility limits of the B- 
and a-phases (Cg a - C a g) . Diffusion in a three-phase system (limited solid 
solubility with intermediate phase formation) produces a concentration profile 
(fig. 1(c)) with discontinuities corresponding to the differences between the 
solubility limits of the 3- and y-phases (Cgy - c yg) and the solubility limits 
of the y- and a-phases (Cy a - C a y) . 

Diffusion in these systems can be described by Fick's second law for each 
phase and by a flux-balance equation for each interface. Fick's second law is 
of the form 


3C 

at 



(i) 


where C is the atomic fraction of B, r is the distance from the center of 
B-rich region, D is the concentration-dependent diffusion coefficient, and 
m = 0, 1, or 2 for planar, cylindrical, or spherical geometries, respectively. 
The interface-flux -balance equation for the aB interface in the two-phase 
system is given by 
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where 5/2 is the location of the aB interface, with superscripts + 
and - designating the a and 3 sides of the interface, respectively. Two 
interface- flux-balance equations are required for the three-phase system. . 
These are of the form 


(C By c y6 ) 


d(5j/2) 

dt 


= D. 
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at the BY interface and 
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at the ay interface, where 5i/2 is the location of the By interface with 
superscripts + and - denoting the y and 3 sides of the interface, 
respectively, and £ 2 /^ is the location of the ay interface with super- 
scripts + and - denoting the a and y sides of the interface, respec- 
tively. The initial and boundary conditions (fig. 2) to be satisfied for all 
three systems are 


C = C 1 


(0 = r < 1/2; t = 0) 


C = C r 


( 1/2 < r ^ L/2; t = 0) 


|£-o 

3r 


(r = 0 and L/2; t ^ 0) 


In this study, C’ and C Q were assumed to be 1 and 0, respectively. 

For two- and three-phase systems, all interface compositions were assumed 
equal to the equilibrium solubility limits of the phases. For a two-phase 
system at aB interface. 


: Ba 

(r = 5V2) 

: a3 

(r = 5 + /2) 
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For a three-phase system at By interface 



and at ay interface, 



(r = q/2) 
(r = Z,\/2) 

(r = Z~/2) 
(r = £ 3 / 2 ) 


Finite-difference analyses and computer programs were developed to solve 
the preceding equations. These solutions are applicable for planar, cylin- 
drical, or spherical geometries with any size for the diffusion zone and any 
continuous variation (within a given phase) of the diffusion coefficient with 
concentration. Special techniques were included in the analyses to account 
for differences in molal volumes, initiation and growth of an intermediate 
phase (three-phase system) , disappearance of a phase (two- and three-phase 
systems) , and the presence of an initial composition profile in the specimen. 

A major improvement in solution accuracy was achieved in the two-phase analysis 
by employing a mass-conservation criterion to establish the location of the 
interface rather than the conventional interface-flux-balance criterion. In the 
three-phase analysis, computation time was minimized without sacrificing solu- 
tion accuracy by treating the three-phase problem as a two-phase problem when 
the thickness of the intermediate phase was less than a small preset value. 

Three computer codes were developed to perform these analyses. Information 
concerning the availability of these codes and the essential features of each 
code, including a discussion of stability and accuracy, is found in reference 13. 

Although the programs are general and can treat volume differences, the 
molal volumes and atomic weights of atoms A and B were assumed to be equal 
for the purposes of this paper. Thus, the concentration of B is the same 
whether expressed as atomic, weight, or volume fraction. 

All three programs treat concentration dependence of the diffusion coeffi- 
cient. Although any continuous variation of D (within a given phase) can be 
treated, only three types of D behavior were considered for this study: 
constant D, log D linearly increasing with C, and log D linearly 
decreasing with C. These cases were selected because of the experimentally 
observed trend for the log D to vary linearly with concentration (ref. 11) . 

The change in D within a phase was assumed to be one of the three types: 
no change, one order of magnitude, and two orders of magnitude. Table 1 shows 
representative diffusion coefficient variations for several single-phase and 
multiphase binary alloy systems. The magnitude and variations of the diffusion 
coefficients chosen for this study are typical of values observed in metallic 
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systems above about two-thirds of the absolute melting point (ref. 14) . The 
alloy phase thicknesses and exposure times were selected to yield concentration 
profiles which are typical of real systems, and to illustrate the effect of 
D variation on the diffusion kinetics. 


RESULTS AND DISCUSSION 
Single-Phase Systems 

Effect of D variation .- Figure 3(a) shows the D variations studied for 
the single-phase system. The diffusion coefficient was assumed to vary lin- 
early, on a log scale, with concentration. Four D variations were considered: 
(1) D B /D A = 0.01, (2) D 8 /^ = 0.1, (3) D B /D A = 1, and (4) D B /D A = 100, 
where D and D are the chemical diffusion coefficients of pure B and 

pure A, respectively. A summary of the conditions for which solutions were 
obtained is given in table 2. In all cases the radius or semithickness (L/2) 
of the couples was 150 pm. 

Figure 3(b) shows the concentration profiles calculated for a cylindrical 
sample with C = 0.111 using D of the form: (1) D B /D A = 100, (2) D®/D A = 1, 
and (3) D B /D A = 0.01 (cases 1, 2, and 3 of table 2) . The exposure time shown 
was selected to give an intermediate stage of diffusion where distinct differ- 
ences in the concentration profiles due to the different D variations were 
readily detectable. In the concentration region above 50 percent B, the con- 
centration gradient dC/dr is largest for the D B /D A = 0.01 curve and smallest 
for the D b /D a = 100 curve. However, below 50 percent B the D B /D A = 100 
variation curve has the largest gradient and the D B /D A = 0.01 curve the 
smallest. The concentration gradient tends to be largest where the diffusion 
coefficients are smallest, as expected. 

Comparison of geometry and D variation effects .- Figure 4(a) shows the 
effect of D variation on the extent of diffusion for cylindrical geometry 
with C = 0.333. The amount of B remaining in the original B-rich region, H, 
is plotted as a function of time for D B /D A = 100, 1, and 0.01 (cases 4, 5, 
and 6 of table 2) . The curve for D B /D A = 0.01 approaches equilibrium value 
fastest because the value of D A is higher than the D A for the other two 
cases. In approaching equilibrium, the D values corresponding to composi- 
tions near pure A are important because of the absence of compositions near 
pure B for this high degree of homogenization. The time required to reach 
H = 0.50 was twice as long for the D B /D A = 100 case as for the D B /D A = 0.01 
case . 


Figure 4 (b) shows the effect of interfacial geometry on the extent of dif- 
fusion for couples with C = 0.333 and D B /D A = 0.10 (cases 7, 8, and 9 of 
table 2) . The effect is as expected with diffusion taking place fastest for 
spherical geometry and slowest for planar geometry. Comparison of the curves 
in figure 4(b) with those in figure 4(a) suggests that the type of D varia- 
tions observed in single-phase systems can be as important in determining the 
kinetics of interdiffusion as is the geometry of the sample. 
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Figure 5 shows the effect of interfacial geometry on the extent of diffu- 
sion for couples with C = 0.500 and 0.010 (cases 10 to 12 and 13 to 15 of 
table 2). For high values of C (fig. 5(a)), diffusion proceeded fastest for 
spherical geometry and slowest for planar geometry. For low values of C 
(fig. 5(b)), just the opposite was true. An explanation of this switchover 
behavior is provided by the calculations presented in figure 6. The ratio s/v 
of interfacial surface area of B-rich material to its volume is shown for each 
geometry. These ratios are expressed in terms of C and L which were held 
constant in figures 5(a) and 5(b). The values of s/v are listed for 
C = 0.500 and 0.010. As can be seen, planar geometry has the highest s/v of 
the three geometries for C = 0.010 and spherical geometry has the highest 
s/v for C = 0.500. 


Two-Phase Systems 

Effect of D variations .- The D variations shown in figure 7 were used 
to determine the effect of a concentration-dependent diffusion coefficient on 
interface motion and concentration profiles in two-phase systems. The D in 
the major-alloy phase a was assumed to vary linearly on a log scale with con- 
centration. An increase in the diffusion coefficient of a in going from the 

interface composition to pure A has been denoted by D a + and a decrease in 

a. 

the diffusion coefficient of a by D +. This notation was adopted because it 
emphasizes the fact that the diffusion coefficients near the interface exert a 
greater influence on the diffusion kinetics than do the coefficients of the 
material away from the interface. 

Two averages of the diffusion coefficient in the a-phase were considered; 
a logarithmic average D a ->, and an integral average <(D a /. These averages were 
calculated as follows: 



where D a g and are the diffusion coefficients corresponding to the compo- 

sition of the interface C a g and to pure A, respectively. 
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The combinations of D and D M for which calculations were performed 
are listed in table 3. The specimen geometry, initial radius (or semithickness) 
of the g-phase region, and average composition of the sample C are also tabu- 
lated for each case considered. The diffusion coefficient in the minor-alloy 
phase 3 was held constant. Also, the solubility limits of the 3- and a-phases, 
*"301 aru ^ c a3' were held constant. 

A plot of the normalized thickness of the 3-phase versus square root of 
the exposure time for C = 0.02 (cases 1, 2, 3, and 4 of table 3) is shown in 
figure 8(a). The fastest rate of interface motion occurred for D a l and the 
slowest for D a + because the diffusion coefficient near the interface is high- 
est for the D a -t' variation and lowest for the D a t variation. The D a -*- vari- 
ation produced interface motion similar to the D a t variation, whereas the 
<D a > produced interface motion more nearly like that of D a +. Thus, the best 
average diffusion coefficient for a D a t variation is a logarithmic average 
(eq. (5)), and for a D a 4' variation, an integral average is best (eq. (6)). 

Figure 8(b) shows concentration profiles at \[t = 300 sec^/^, calculated 
for cases 1, 2, and 4 of table 3. Saturation of the 3-phase occurs early 
because the 3"Ph ase I s thin relative to the a-phase and the diffusion coeffi- 
cients are approximately the same in each phase. After the 3-phase saturates, 
it will decrease in thickness at a rate controlled by diffusion in the a-phase. 
The concentration gradient in the a-phase near the interface is steepest for 
D a t because of the low D a g. The curves for D 01 -)- and D a t crossover away 
from the interface because, as the concentration tends to pure A, the diffusion 
coefficients for D a + are higher than those for D a l. 

The extent of interdiffusion for a cylindrical couple with C = 0.20 is 
shown in figure 9 where £/£ and H (normalized amount of B in the 3-phase) 
are plotted as a function of \ft. Both plots show essentially the same behavior 
with equilibrium being established at a £/£ value of 0.598 and an H value 
of 0.304. Both a- and 3~phases are present in the sample at equilibrium 
because C = 0.20 lies in the a + 3 region. 

Comparison of geometry and D variation effects .- A comparison of the 
effects due to differences in geometry and in D variations for a two-phase 
system is shown in figure 10. The results indicate that typical D a varia- 
tions can have a greater influence on the kinetics of interdiffusion than the 
geometry of the interface. Comparison of the results of figure 10 for a two- 
phase sample with the results of figure 4 for a single-phase sample reveals 
that D variations in two-phase systems can be as important as D variations 
in single-phase systems . 


Three-Phase Systems 

The type of diffusion behavior observed in three-phase systems is depen- 
dent on solubilities of each phase, rates of diffusion in each phase, and 
initial thicknesses of the a- and 3-phases, and therefore the average composi- 
tion of the diffusion couple C. If C lies in one of the single-phase 
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regions, only that phase will remain at equilibrium. If C lies in the a + Y 
region, the equilibrium volume fractions of y- and a-phases are given by 



Y C - C 
ya ay 

and 

V* = 1 - V* 

a y 

In this study, C was chosen to be in the a region or the a + y region. 

Concentration profiles .- The principal effects of the solubility limits 
and C are, to a large extent, predictable from the equilibrium phase diagram. 
However, the effects of unequal rates of diffusion in each phase and variations 
in the diffusion coefficient within a given phase are not discernible from the 
phase diagram. For this reason, the primary emphasis of this study was placed 
on characterizing the effect of the different diffusion coefficients and their 
variations on the dissolution characteristics. The D variations investigated 
in this study are illustrated in figure 11. The different D combinations con- 
sidered are listed in table 4 along with C and the solubility limits. These 
cases were selected because they represent typical parameters for real systems. 
The designations for D variations (D D f, D +, and <D>) are the same as 

the ones discussed for two-phase systems. 

The concentration profiles calculated for a cylindrical couple with 
C = 0.02 (case 1 of table 4) are shown in figure 12. The three profiles 
illustrate the major concentration and phase changes occurring during the 
initial and intermediate stages of diffusion. Because the B-phase is thin 
relative to the a-phase, diffusion proceeds to the center of the B-phase early 
in the diffusion process resulting in a reduction in the concentration gradient 
in the B _ phase. The intermediate y-phase does not significantly increase in 
thickness until the concentration gradient in the B~phase approaches zero. 

From this time on, the y-phase grows at the expense of the B~phase until the 
B-phase is totally consumed. The y-phase will then decrease in thickness and 
disappear with only the a-phase remaining at complete homogenization. The 
sequence of events is controlled by solubilities, initial volume fractions, and 
diffusion coefficients. The effect of each of these variables on the diffusion 
process will be discussed in the following sections with primary emphasis 
placed on the effect of diffusion coefficient variations on the motion of 
phase boundaries . 

Effect of relative magnitude and variation of D in the major-alloy 

a a a / a \ B 

phase .- Solutions were determined for D + , Df, D -s-, and \D / assuming D 
and to be constant and equal to D a ->- (cases 1, 2, 3, and 4 of table 4) . 

The results for the normalized thickness of the B-phase are shown in figure 13 
as a function of \[t. The curves in this figure have the same general behavior 
as those presented in figure 8(a) (same C, geometry, D variations, and 
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terminal phase solubilities) for the two-phase sample. The principal differ- 
ence in the two sets of curves is that, in the three-phase system, the 6-phase 
disappears in a shorter: time and the spread between the D variation curves 
is not as large. The faster rate of loss of 6 in the three-phase system is 
due to its low solubility gap (<$gy = 0.85 - 0.51 = 0.34) when compared to the 
two-phase value (6^ = 0.85 - 0.15 = 0.70), requiring less material transport 
for the 6 to y transformation with respect to 6 to a transformation. 

The spread between the .D variation curves in a three-phase system is smaller 
because the intermediate y-phase damps the effect of the D a variation. This 
damping effect coupled with lower solubility limits makes the D a variations 
in three-phase systems less important than in two-phase systems. 

The sensitivity of the diffusion process to large differences in the value 
of D a (with respect to and dY) is illustrated in figure 14. The con- 

ditions for these calculations are summarized in table 4 (cases 5, 6, and 7). 
When D a was high ((D a ->-) x 100), the thickness of the 6-phase decreased rapidly 
until the a-phase saturated (fig. 14(a)). The concentration profiles shown in 
figure 14(b) indicate that at \ft = 50 sec'''^, the a-phase is saturated for the 
(D a ->-) x 100 case. Once the a-phase saturates, the kinetics are primarily con- 
trolled by diffusion in the y-phase because the solubility of 6~phase is only 
1 percent. 

Effect of relative magnitude of D in minor-alloy phase .- Because D 
variations in the minor-alloy phase are less important than in the major-alloy 
phase and would, in the worst case, have only an initial transient effect, only 
order of magnitude changes in the diffusion coefficient of the minor-alloy 
phase were studied (cases 1, 8, and 9 of table 4) . The normalized thicknesses 
of the 6~ and y-phases are shown in figure 15 as a function of \ft. The change 
in thickness of the y-phase, (£ 2 - £^)/£, was not significantly affected by the 
6 6 

value of D. The value of D has only an initial transient effect on the 
thickness of the 6~phase, (^y/£) . Once the 6-phase saturates, the value of 
is unimportant because and D a then control the diffusion process. 


Effect of relative magnitude of D in the intermediate phase .- Because 
intermediate phases in binary alloy systems generally have small solubility 
ranges (less than 0.05, atomic fraction B) , concentration dependence of dif- 
fusion coefficients in these phases is not significant. Therefore, in studying 
the effect of the diffusion coefficient of the intermediate phase dY on the 
diffusion kinetics, only order of magnitude changes in dY were considered. 
Figure 16 shows the results for cases 1, 10, and 11 of table 4. When dY was 
low ( (dY->) / loo) , the y-phase did not grow in thickness, and the rate of loss 
of 6-phase was less than for D^-v. Under these conditions, computation time is 
minimized (without sacrificing solution accuracy) by analytically treating the 
three-phase problem as a two-phase problem. When was high ((DY->-) x 100), 

the y-phase grew rapidly at the expense of both a and 6; this resulted 
in a rapid decrease in the 6-phase thickness. Concentration profiles for 
\ft = 50 sec 1 ' are shown in figure 16(b) for all three values. The pro- 

file for (D^-*-) x 100 shows that the y-phase grew into the initial 6-phase 
region much faster than into the a-phase because of the early saturation of the 
6-phase. The y-phase will increase in thickness until 6 disappears and then 
decrease in thickness until it disappears leaving only a-phase in the specimen. 
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Figure 17 shows the results for C = 0.40 (cases 5, 12, and 13 of 
table 4) . Although C and solubilities are different from those used to cal- 
culate the curves shown in figure 16, the same general behavior is observed. 

If dY is less than or equal to D a and D^, the kinetics are primarily con- 
trolled by diffusion in the a-phase. However, if is larger than D a 

and D p , the kinetics are dominated in the early stages by the growth of the 
y-phase. Therefore, intermediate phases in three-phase , systems are unlikely to 
act as suitable diffusion barriers between the terminal phases. 


CONCLUDING REMARKS 

Diffusion calculations were performed to establish the conditions under 
which concentration dependence of the diffusion coefficient was important in 
single-, two-, and three-phase binary alloy systems. Finite-difference solu- 
tions for each type of system were obtained using diffusion coefficient varia- 
tions typical of those observed in real alloy systems. Solutions were also 
obtained using average diffusion coefficients and used as reference in assessing 
the effects of diffusion coefficient variations. Calculations were performed 
for planar, cylindrical, and spherical geometries in order to compare the effect 
of diffusion coefficient variations with the effect of interface geometries. 
Based on the results obtained, the following conclusions were drawn. 


Single-Phase Systems 

1. Concentration dependence of the diffusion coefficient in single-phase 
systems was found to affect the kinetics of diffusion as strongly as the inter- 
facial geometry of the diffusion couple. 

2. The influence of geometry and average composition of the couple 

(C, atomic fraction B) on diffusion kinetics is related to the ratio of inter- 
facial surface area of B-rich material to its volume. 


Two-Phase Systems 

1. Concentration dependence of the diffusion coefficient in the major- 
alloy phase of a two-phase system affects the kinetics of diffusion as strongly 
as the interfacial geometry of the diffusion couple. 

2. The diffusion coefficients of the compositions close to the solubility 
limits exert a greater influence on the diffusion kinetics than do the coeffi- 
cients of the compositions further away from the interface. 

3. The integral average gives a better representation of the major-alloy 
phase diffusion coefficient D 01 when D a decreases in going from the inter- 

facial composition toward pure material. The logarithmic average gives a 
better representation of D 3 * 05 variation when D a increases in going from the 

interfacial composition toward pure material. 
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Three-Phase Systems 


1. Concentration dependence of the diffusion coefficient is significant 
only for the major-alloy phase. The damping effect of the intermediate phase 
makes the effect of diffusion coefficient variations on minor-phase thickness 
less important with respect to the two-phase system for comparable conditions. 

2. Two orders of magnitude change in the diffusion coefficient of the 
minor-alloy phase for a sample with C = 0.02 had only an initial transient 
effect on the kinetics of diffusion. 

3. The diffusion coefficient of the intermediate phase has significant 

effect on the diffusion kinetics only when is greater than the diffusion 

coefficients of the terminal phases. 

4. When is small relative to the diffusion coefficients of the termi- 

nal phases, the intermediate phase remains extremely thin during the initial 
stages of diffusion. Under these conditions, computation time is minimized 
(without sacrificing solution accuracy) by analytically treating the three-phase 
problem as a two-phase problem. 

5. Intermediate phases in three-phase systems are unlikely to act as 
suitable diffusion barriers between the terminal phases. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
August 25, 1978 


13 



REFERENCES 


1. Tanzilli, R. A.; and Heckel, R. W. : Numerical Solutions to the Finite, 

Diffusion-Controlled, Two-Phase, Moving-Interface Problem (With Planar, 
Cylindrical, and Spherical Interfaces) . Trans. Met. Soc. AIME, vol. 242, 
no. 11, Nov. 1968, pp. 2313-2321. 

2. Tanzilli, R. A.; and Heckel, R. W. : An Analysis of Interdiffusion in 

Finite-Geometry, Two-Phase Diffusion Couples in the Ni-W and Ag-Cu 
Systems. Metall. Trans., vol. 2, no. 7, July 1971, pp. 1779-1784. 

3. Lanam, R. D.; and Heckel, R. W. : A Study of the Effect of an Intermediate 

Phase on the Dissolution and Homogenization Characteristics of Binary 
Alloys. Metall. Trans., vol. 2, no. 8, Aug. 1971, pp. 2255-2266. 

4. Unnam, J.; and Houska, C. R. : Iterative and Finite Difference Computer 

Programs for Diffusion in Single and Two Phase Systems. Grant 

No. NSF-GH-30180A #1, Virginia Polytech. Inst. & State Univ. , Jan. 1976. 

(Available from NTIS as PB-256 282.) 

5. Houska, C. R. ; and Unnam, J. : Iterative Approaches Describing Atomic Dif- 

fusion in Finite Single- and Two-Phase Systems. J. Appl. Phys., vol. 47, 
no. 10, Oct. 1976, pp. 4325-4335. 

6. Herring, Harvey W. ; and Tenney, Darrel R. : Diffusion in a Unidirectional 

Filament Reinforced Metal Composite. Metall. Trans., vol. 4, no. 2, 

Feb. 1973, pp. 437-441. 

7. Tenney, Darrel R. : Diffusion Analysis for Two-Phase Metal-Matrix Composite. 

NASA TN D-8299 , 1976. 

8. Tenney, Darrel R. ; and Talty, Patrick K. : X-Ray Diffraction Investigation 

of Bimetallic Diffusion Zones in the Cu-Pd System. Metall. Trans., 
vol. 5, no. 1, Jan. 1974, pp. 241-247. 

9. Unnam, J.; Carpenter, J. A.; and Houska, C. R. : X-Ray Diffraction Approach 

to Grain Boundary and Volume Diffusion. J. Appl. Phys., vol. 44, no. 5, 
May 1973, pp. 1957-1967. 

10. Unnam, J. ; and Houska, C. R. : An X-Ray Study of Diffusion in the Cu-Ag 

System. J. Appl. Phys., vol. 47, no. 10, Oct. 1976, pp. 4336-4342. 

11. Smithells, Colin J. : Metals Reference Book - Volume II. Fourth ed. 

Plenum Press, 1967. 

12. Svechnikov, V. N. , ed. : Diffusion Processes in Metals. TT 70-50030, 

Nat. Bur. Stand, and Nat. Sci. Found., c.1970. 


14 



13. Tenney, Darrel R. ; and Unnam, Jalaiah: Numerical Analyses for Treating 

Diffusion in Single-, Two-, and Three-Phase Binary Alloy Systems. NASA 
TM-78636, 1978. 

14. Seith, Wolfgang; and Heumann, Theodor: Diffusion of Metals: Exchange 

Reactions. Second ed., AEC-tr-4506, U.S. At. Energy Comm., 1955. 


15 



TABLE 1.- DIFFUSION COEFFICIENT VARIATIONS FOR SEVERAL 
SINGLE-PHASE AND MULTIPHASE BINARY ALLOY SYSTEMS 


System, 

A-B 

Temp . , 
K 

Composition 
range , 
atomic 
fraction B 

Cu-Ni 

1173 

0 - 1.00 

Ni-Pd 

1292 

0 - 1.00 

Cu-Pd 

1173 

0 - 1.00 

Nb-W 

2373 

0 - 1.00 

Ag-Cu 

1023 

0-0.12 Cu 

Mo-Ti 

1123 

0.10-0.90 Ti 

Ag-Zn 

883 

0.40-0.55 Zn 

Fe-Al 

1123 

0-0.18 A1 


Max. D 


Min. D 


No. of 
phases 


Reference 


12 (p. 52) 

8 

11 (p. 679) 
10 

12 (p. 39) 

11 (p. 674) 
11 (p. 675) 


















No 


Geometry 




C 


( 1 / 2 ) , ym 


1 

Cylindrical 

100 

0.111 

50 

2 

Cylindrical 

1 

.111 

50 

3 

Cylindrical 

.01 

.111 

50 

4 

Cylindrical 

100 

.333 

86.6 

5 

Cylindrical 

1 

.333 

86 . 6 

6 

Cylindrical 

.01 

.333 

86.6 

7 

Cylindrical 

.10 

.333 

86.6 

8 

Planar 

.10 

.333 

50 

9 

Spherical 

.10 

.333 

104 

10 

Planar 

1 

.500 

75 

11 

Cylindrical 

1 

.500 

106.1 

12 

Spherical 

1 

.500 

119.1 

13 

Planar 

1 

.010 

1.5 

14 

Cylindrical 

1 

.010 

15 

15 

Spherical 

1 

.010 

32.3 










TABLE 3.- CASES CONSIDERED FC 
fc ft = 0.85, Co = 0.15, 


Geometry 

D a 


Cylindrical 

D a i 

D^-> 

Cylindrical 

D a i 

D^-> 

Cylindrical 

D a -> 

D^-s- 

Cylindrical 

<D a > 

D^-> 

Cylindrical 

D a i 

D&->- 

Cylindrical 

D a + 

D&- 

Cylindrical 

D 01 -)- 

D^-> 

Cylindrical 

<D a > 

D$->- 

Planar 

rv ! 

D U -> 

D^+ 

Spherical 

D a + 

A 


L TWO-PHASE SYSTEM 
L/2 = 140 ym] 


C 

ii/2 ) , ym 

0.02 

20.0 

.02 

20.0 

.02 

20.0 

.02 

20.0 

.20 

62.6 

.20 

62.6 

.20 

62.6 

.20 

62.6 

.20 

28.0 

.20 

81.9 













TABLE 4.- CASES CONSIDERED FOR THREE-PHASE SYSTEM 
[cylindrical geometry, c yb = °- 51 ' V - °- 49 ] 


No. 

D a 

D^ 

d y 

C $Y 

c 

u ay 

C 

( V2) , ym 

(L/2) , 

1 

D a + 

D^->- 

d y -* 


0.85 

0.15 

0.02 

20 

140 

2 

D 0t 4' 

D^ 

d y -> 


.85 

.15 

.02 

20 

140 

3 

D a t 

D^-J- 

d y -> 


.85 

.15 

.02 

20 

140 

4 

<D a > 

D^->- 

d y + 


.85 

.15 

.02 

20 

140 

5 

D a -v 

D^ 

d y + 


.99 

.20 

.40 

50 

79 

6 

(D a ->) x 100 

D^ 

d y + 


.99 

.20 

.40 

50 

79 

7 

(D a +)/100 

D®-* 

d y + 


.99 

.20 

.40 

50 

79 

8 

D a + 

r\j 

(D^->) x 100 

d y ->- 


.85 

.15 

.02 

20 

140 


D a ->- 

(D P ->-)/100 

d y + 


.85 

.15 

.02 

20 

140 

9 

D a + 

D^-> 

(D Y ->-) x 

100 

.85 

.15 

.02 

20 

140 

B 

D 01 -, 

D^-> 

(d y ->-) /100 

.85 

.15 

.02 

20 

140 

12 

D a -> 

D^ 

(D Y ->-) X 

100 

.99 

.20 

.40 

50 

79 

13 

D a ^ 

D^-»- 

(D Y -0/100 

.99 

.20 

.40 

50 

79 



















(a) Single phase - complete solid solubility. 



(b) Two phase - limited solid solubility. 



(c) Three phase - intermediate-phase formation. 


Figure 1.- Phase diagrams and examples of concentration profiles produced by 
diffusion between initially pure A and pure B at temperature 0 D in 
single-, two-, and three-phase binary alloy systems. 
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Planar 


Cylindrical 


Figure 2 



C = (£/L ) 2 


- Average composition, initial conditions, and zero-flux (J = 0) 
boundary conditions for different geometries. 
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Atomic percent B 


Atomic percent B 

(a) D variations. 


Initial Interface 


(b) Concentration profiles. 

Figure 3.- Example of D variations in single-phase systems and associated 
concentration profiles in cylindrical couple with C = 0.111 at inter- 
mediate stage of diffusion (5 hr) . 













Figure 5.- Average composition and geometry effects on 
degree of homogenization. 
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s _ lxl 
v f/2 x 1 x 1 

C = (A/L) 

_s_ = 2 

v CL 


C (s/v) X L/2 = 


7T a X 1 

7T (C/2) 2 X 1 

= (*/L) 2 

= 4 

CC)V2 L 


2.83 

20.0 



V 

4/3 7 T (St/2) 3 1 

C = 

(J/L) 3 

s 

6 

V 

(C) 1 / 3 L 

(s/v) 

X L/2 = 37(C) 1 / 3 


3.78 


13.9 
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(a) Change in thickness of 8-phase with diffusion time. 



Vt. sec 1/2 

(b) Change in amount of B in 8-phase with diffusion time. 


Figure 9.- Effect of D a variations in cylindrical couple with C = 
Diffusion conditions correspond to cases 5, 6, 7, and 8 of table 
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Planar 

Cylindrical D c 
Spherical 


Equilibrium Value 


Figure 10.- Geometry and D variation effects on degree 
of homogenization. Diffusion conditions for D varia 
tion curves correspond to cases 5, 6, and 7 of table 3 
Conditions for different geometry curves correspond to 
cases 7, 9, and 10 of table 3. 





Atomic fraction B 

Figure 11.- Diffusion coefficient variation in a-, and yphases 

of three-phase system. 
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Atomic percent B 



different times in cylindrical 
conditions correspond to case 1 








0 50 100 150 200 250 



(a) Change in thickness of y- and B-phases 
with diffusion time. 



0 20 40 60 80 


r. Mm 

(b) Concentration profiles at \[t = 50 sec^/^ . 

Figure 14.- Effect of D a variations in cylindrical couple with 
C = 0.40 on y- and B-phase thicknesses and concentration 
profiles. Diffusion conditions correspond to cases 5, 6, 
and 7 of table 4 . 











Atomic percent B 


(a) Change in thickness of y- and B-phases with diffusion. 



r , Mm 

(b) Concentration profiles at \[t = 50 sec 1 / 2 . 

Figure 16.- Effect of large changes in for cylindrical couple with 

C = 0.02 on Y- and B.-phase thicknesses and concentration profiles. 
Diffusion conditions correspond to cases 1, 10, and 11 of table 4. 
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